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The identity and characteristics of the lifetime limiting de¬ 
fects in n-type 4H-SiC epitaxial layers are of particular cur¬ 
rent interest, due to the suitability of this material for high- 
power, solid-state switching devices. Much work has been 
done in the past decade to identify the spectral signature and 
the local atomic structure of the defect that controls the life¬ 
time in this material. Until recently, it was concluded that 
two dominant electron traps, Z 1/2 and EH 6/7 , both associated 
with a carbon vacancy, were controlling the lifetime. DLTS 
and optical studies show, however, that EH 6/7 is not an impor¬ 


tant minority carrier trap, due to its small capture cross sec¬ 
tion for holes. Thus, Z]/ 2 acts alone as the lifetime limiting 
defect. The details of recombination at this trap have been 
studied through the injection level and temperature depend¬ 
ence of the carrier lifetime and through carrier dynamics 
simulations that take into account the complex nature of the 
defect. The simulations shed light on the processes that domi¬ 
nate the recombination at Z^over a large range of injection 
level and temperature. 
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1 Introduction Wide bandgap semiconductors are of 
great current interest, as their material properties - high 
thermal conductivity, high breakdown voltage, high ther¬ 
mal stability and chemical inertness - make them ideal for 
applications that demand high power and/or high fre¬ 
quency operation that is well beyond the limitations of 
Si-based technology. Consequently, there is much interest 
in these systems for high power solid-state switching 
and communications applications. As a result, much re¬ 
search has been carried out in the direct-gap GaN and 
GaN/AlGaN systems for high power, high frequency mi¬ 
crowave transistors and in SiC (indirect band gap) for high 
power bipolar switching devices. 

SiC grows in numerous possible stacking orders 
(polytypes), each with a distinct bandgap. The most com¬ 
mon polytypes are 3C (cubic, zincblende structure, 
E g = 2.36 eV), 6H and 4H (hexagonal, wurtzite structure, 
7i g = 3.0eV and 3.23 eV, respectively). For applications 
involving high-voltage bipolar switching devices, it is clear 
from a comparison of typical material parameters for GaN 
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and SiC to those of Si, shown in Table 1, that the combina¬ 
tion of high thermal conductivity and high breakdown field 
provides SiC with a distinct advantage over both GaN and 
Si. The 4H polytype of SiC is favored over 6H in these 
high power vertical devices because of its much higher 
mobility along the c-axis, which is the growth direction in 
most epitaxial layers. 

High power switching devices, such as p-i-n diodes 
and insulated-gate bipolar transistors, require a thick, low- 
doped drift layer to provide the low conductivity necessary 
to hold off a high voltage across the device when the 
switch is in the off-state. This is facilitated by the high 
breakdown field of SiC. With the switch in the on-state, 
however, the same n drift region must have a high con¬ 
ductivity, in order to assure a low forward voltage drop and 
low power dissipation. This duality in the conductivity of 
the drift layer is made possible by the process of conduc¬ 
tivity modulation, whereby a high concentration of minor¬ 
ity carriers (and majority carriers) are injected into the drift 
region under high forward bias, resulting in a large en- 
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Table 1 Comparison of material properties of Si, SiC and GaN. 
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3.39 
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3.3 
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hancement in the conductivity of the layer. The degree to 
which conductivity modulation is effective in providing the 
necessary high carrier concentration to the drift layer is di¬ 
rectly related to the lifetime of the minority carriers within 
the region: a short carrier lifetime results in limited con¬ 
ductivity modulation and a large forward voltage drop. 
Consequently, understanding the limitations on the carrier 
lifetime has become a topic of intense interest. The carrier 
lifetime in the drift region affects not only the conductivity 
of the layer under forward bias, but also the switching 
speed of the device: Longer lifetimes afford better conduc¬ 
tivity modulation but also slower switching speeds. Conse¬ 
quently, device design requires a compromise between 
these two parameters, so it is necessary to be able to tailor 
the carrier lifetime for a particular device. This kind of life¬ 
time control can be achieved if the minority carrier lifetime 
in the as-grown drift layer can be made longer than re¬ 
quired and if a post-growth process, such as electron irra¬ 
diation, is developed to reduce the lifetime in a controlled 
manner to the required level [1], Both conductivity modu¬ 
lation and lifetime control require long carrier lifetimes. 

The primary mechanism responsible for reducing the 
carrier lifetime is the trapping of free carriers at material 
defects. Consequently, it has become important to identify 
and to characterize the specific defects that dominate the 
trapping of minority carriers and thus control the carrier 
lifetime within the drift region. In Section 2 we review 
briefly the evolution, over the last 10-15 years, of a) the 
microscopic structures that have been associated with the 
dominant defects, and b) the correlation between the car¬ 
rier lifetime and specific defects. The progress toward 
identifying the responsible recombination centers led to the 
conclusion that the lifetime limiting defects were the two 
dominant electron traps observed in deep level transient 
spectroscopy (DLTS) measurements of as-grown epilayers, 
labeled as Z 1/2 (or Z1/Z2) and EH 6/7 (or EH6/EH7). How¬ 
ever, evidence presented in Section 3 indicates that it is 
only Z, /2 that has a significant effect upon the carrier life¬ 
time. The processes controlling the carrier dynamics at this 
defect are considered in detail in Section 4, both through 
experiment and through carrier dynamics simulations, in 
order to provide a deeper understanding of recombination 
at this defect and of those parameters that significantly af¬ 
fect the lifetime. The discussion is summarized in Sec¬ 
tion 5. 



0 200 400 600 800 


Temperature (K) 

Figure 1 DLTS (a) and MCTS (b) spectra of typical n-type 
4H-SiC epitaxial layers. (Adapted from Ref. [2].) 

The DLTS spectrum of a typical as-grown, n-type epi¬ 
taxial layer of 4H-SiC, similar to that shown in Fig. 1(a) 
[2], is dominated by the two electron traps, generally re¬ 
ferred to as Z 1/2 (~Zi c - 0.65 eV) and EH 6/7 (~E C - 1.65 eV). 
Each of these designations refers to two similar but distin¬ 
guishable defects, Z ls Z 2 and EH 6 , EH 7 - each pair thought 
to be the same defect on cubic and hexagonal lattice 
sites [3], Other traps are often observed, such as Ti 
( E c - 0.17 eV), but usually at lower concentrations. Minor¬ 
ity carrier transient spectroscopy (MCTS) measurements in 
as-grown material, similar to that shown in Fig. 1(b) [2], 
are generally dominated by two persistent boron-related 
hole traps: the shallow boron acceptor (£ v + 0.27 eV) and 
the D-center(E v + 0.65 eV), a boron-related complex [4, 5], 
DLTS, in conjunction with carrier lifetime measurements, 
has been central in the effort to identify the defects that 
control the carrier lifetime. Much of this work has been fo¬ 
cused on either a) identifying the atomic structure of im¬ 
portant traps, by correlating DLTS features with varying 
growth or post-growth processing conditions, or b) identi¬ 
fying the traps that limit the carrier lifetime, through corre¬ 
lations between DLTS measurements and carrier lifetime 
measurements. Initially, these two approaches were taken 
almost independently, with significant overlap occurring 
only later. Consequently, the discussion of this work will 
treat each approach separately, beginning with the former. 

2 Defects in 4H-SiC epilayers 

2.1 Defect structure The application of DLTS 
measurements to the study of deep traps in 4H-SiC was 
first carried out on bulk n-type samples grown by sublima¬ 
tion vapor transport by Uddin et al. [6], The spectra were 
dominated by a broad peak due to an electron trap near 
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E c - 0.72 eV, which was thought to be associated with a 
defect at both cubic and hexagonal sites. DLTS measure¬ 
ments on epitaxial layers of n-type 4H-SiC were initially 
investigated by Kimoto et al. [7], who observed a similar, 
dominant electron trap, ~0.61-0.68 eV below the conduc¬ 
tion band. This defect was noted to be similar to the Z,/Z 2 
center reported in 6H-SiC by Pensl and Choyke [8] and by 
Zhang [9, 10], This designation was adopted in 4H-SiC as 
well. Doyle et al. [11] reported a similar dominant trap in 
as-grown 4H-SiC that evolved into two closely-spaced lev¬ 
els after electron irradiation. A similar broadening of the Z, 
DLTS spectrum after He + ion implantation, due to the gen¬ 
eration of a second, similar defect, was observed by Dali- 
bor et al. [12, 13], who found that the Z, component be¬ 
haved as an acceptor and was found to be thermally stable 
above 2000 °C, while the second component annealed out 
near 1000 °C. This was interpreted as evidence that the lat¬ 
ter component was a single intrinsic defect. Using the pre¬ 
viously published association of the Z-center with the D1 
photo luminescence (PL) band [14] in 6H-SiC, the authors 
concluded from low temperature PL and DLTS studies that 
the Z 1/2 center was consistent with a nearest-neighbor 
V c -V si divacancy. However, the Dl-Z 1/2 association in 
4H-SiC was later shown to be erroneous [15], The EH 6/7 
defect was introduced by Hemmingsson et al. [16], appear¬ 
ing as a broad peak (reflecting two defects, EH 6 and EH 7 ) 
in the DLTS spectra of as-grown and electron-irradiated 
samples, with an acceptor level at E c - 1.65 eV and a large 
electron capture cross-section (<r n > 5 x 1CL 15 cm 2 ). 

The broad DLTS signature for Z 1/2 was shown by Hem¬ 
mingsson et al. [3, 17] to result from the fact that the defect 
is a negative-U center and thus able to trap two electrons. 
This results in three possible charge states (empty, singly- 
and doubly-occupied). Consequently, there are two transi¬ 
tions within the gap; a (0/+) donor state and a (-/0) accep¬ 
tor state that lies deeper in the band gap, due to large lattice 
relaxation around the negatively charged center. As a result, 
during a standard DLTS temperature scan thermal ioniza¬ 
tion of an electron from the doubly-occupied state led 
to the immediate ionization of the resulting, shallower “0” 
state, so that an effective -/+ transition was observed. 
Hence, only the single, broad Z 1/2 signature reflecting the 
-/+ transition was seen in DLTS, as shown [3] in Fig. 2(a). 
The signature of the “0” state was not observed. By ini¬ 
tially emptying the ” state with a light pulse and using 
very short filling pulse widths to minimize any initial cap¬ 
ture of carriers back into the ” state, the DLTS signature 
of the 0/+ transition for both Z, and Z 2 were revealed, as 
shown in Fig. 2(b). The configuration coordinate diagram 
of the three Z I/2 charge states in the inset reflected the large 
lattice relaxation of the ” state. 

While several defects were reported from DLTS stud¬ 
ies, there was little progress in identifying the chemical na¬ 
ture of these traps until systematic studies were carried out 
to determine the effect of varying growth conditions and 
post-growth processing on defect concentrations. In the 
case of Z 1/2 , this led to a long series of contradictory identi- 



Figure 2 DLTS specta (a) of the Z 1/2 (—/+) transition, and (b) of 
the Z 1/2 (0/+) transition, obtained by pre-irradiating the sample 
and using short filling pulses. The configuration coordinate dia¬ 
gram in the inset reflects the negative-U property of the defect. 
(Adapted from Ref. [3].) 

fications. By investigating the effect of a varying C/Si ratio 
on the formation of Z x defects during CVD growth, Ki¬ 
moto et al. [18] observed a decrease in the defect concen¬ 
tration for C-rich growth conditions, concluding that Z, 
was an intrinsic defect, most probably a Si antisite (Si c ) or 
C vacancy (V c ). On the other hand, Kawasuso et al. [19] 
compared the concentrations of Z 1/2 defects, measured by 
DLTS, and vacancies, detected by positron annihilation 
studies, as a function of thermal annealing temperature. 
The similarity in their annealing behavior led to the con¬ 
clusion that the defect was related to Si vacancies (V si ). 
Another alternative origin for the Z 1/2 center was suggested 
from the work of Pintilie et al. [20], where the defect con¬ 
centration was positively correlated with the nitrogen dop¬ 
ing level and the C/Si ratio. However, the Z 1/2 concentra¬ 
tion was observed to increase with more C-rich growth 
conditions, which is opposite to previous [18] and later re¬ 
ports. As a result of these investigations, the Z 1/2 center was 
assigned to a substitutional N - interstitial C complex, 
N-C;. This view was supported by theoretical calculations 
[21] suggesting that the Z 1/2 center was consistent with an 
interstitial N - interstitial C complex, Nj-C;. However, the 
involvement of N in the composition of the Z 1/2 defect was 
later refuted [22-24], 

A detailed study of the effect of varying growth condi¬ 
tions (growth temperature, C/Si ratio and sample position 
relative to the gas flow) on the incorporation of defects 
was carried out by Zhang et al. [25] for the Ti and B impu¬ 
rities, the HS1 hole trap and the Z 1/2 and EH 6/7 electron 
traps. The latter two defects exhibited a very similar de¬ 
pendence on growth temperature (with a formation energy 
near 4 eV), similar variations in concentration along the di- 
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C/Si ratio 

Figure 3 Similar dependences of the Z 1/2 and EH 6/7 concentra¬ 
tions on C/Si ratio. (Adapted from Ref. [27].) 

rection of gas flow and a relatively weak dependence on 
C/Si ratio over the range that was studied. From these re¬ 
sults it was concluded that the Z 1/2 center was consistent 
with a divacancy, an antisite pair or an antisite-vacancy 
pair, while the EH 6/7 defect was tentatively associated with 
the carbon vacancy. Kimoto et al. [26] found a very differ¬ 
ent dependence on C/Si ratio for EH 6/7 and Z 1/2 , which was 
consistent with their earlier work: The concentrations of 
both defects decreased significantly as the growth became 
more C-rich. This dependence upon C/Si ratio was later 
shown [27, 28] to be essentially the same for both defects, 
as shown [27] in Fig. 3. The two defects were also found 
[29, 30] to exhibit a very similar dependence upon thermal 
anneal temperature (Fig. 4) - both were found stable to at 
least 1600 °C and decreased by an order of magnitude at 
higher temperatures. All of these similarities suggested that 
the two defects might both be related to the same center, 
such as a V c . 

Storasta et al. [23] investigated the defect formation 
and annealing properties of samples exposed to low energy 
electron irradiation that was below the threshold for Si 



Figure 4 Similar dependences of the Z 1/2 and EH 6/7 concentra¬ 
tions on anneal temperature. (Adapted from Ref. [29].) 


atom displacement. Consequently, only intrinsic defects re¬ 
lated to the initial displacement of carbon atoms from their 
lattice sites were created. The irradiation-induced defects 
that resulted were Z 1/2 and EH 6/7 , as well as the EH1 and 
EFI3 electron traps and the HS2 hole trap. Because only 
carbon atoms were displaced, it was concluded that the ap¬ 
pearance of these defects was found consistent with struc¬ 
tures involving V c , Q C si and Si c , while divacancies, di¬ 
interstitials and impurity complexes (including N-defect 
complexes) were ruled out. These results put strict limita¬ 
tions on the microscopic structures that were possible for 
these defects. Danno et al. [31] suggested that the charac¬ 
teristics of the EH 6/7 defect were consistent with the accep¬ 
tor state of the V c , and that the isolated V c may be associ¬ 
ated with the PI hole trap in p-type material. A further 
study by Litton et al. [32] of the effect of C/Si ratio on Z 1/2 
concentration concluded that the defect was either V c or 
Si c -related. 

Overall, the evidence described above, developed over 
a substantial period of time, indicated that both Z 1/2 and 
EFI 6/7 appear to behave similarly with varying growth con¬ 
ditions (C/Si ratio, growth temperature) or post-growth 
processing (thermal anneal, electron irradiation), suggest¬ 
ing that both centers may be related to the same intrinsic 
defect. This picture has been confirmed recently [33], veri¬ 
fying many of these observations within the same group of 
samples and showing (Fig. 5) that, for a large number of 
as-grown, electron-irradiated and irradiated/annealed sam¬ 
ples, there is a one-to-one correlation between the concen¬ 
trations of Z 1/2 and EFI 6/7 . As both defects have also been 
shown to be consistent with the involvement of a carbon 
vacancy, the current view is that both of these defects are 
V c -related centers. 

It was also noted [33] that while the energy levels 
within the band gap for these two defects were shown to be 
consistent with the energy level positions of different 
charge states of the isolated V c , determined from theoreti¬ 
cal calculations [34, 35] and by photo-EPR measurements 
[36], the involvement of a V c -related complex could not be 



Figure 5 Correlation between Z, /2 and EH 6/7 concentrations for a 
large range of 4H-SiC epilayers. Reprinted with permission from 
Ref. [33], © 2006, American Institute of Physics. 
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ruled out. The picture becomes somewhat more complex 
when the correlation between calculated defect properties 
and experimental results includes, in addition to energy 
levels, the donor or acceptor nature of the defect states and 
whether the defect forms a —U center. With these consid¬ 
erations in mind, Eberlein et al. [37] concluded that while 
both Z 1/2 and EH 6/7 could be consistent with different V c 
charge states, other possibilities exist, such as the assign¬ 
ment of Z 1/2 as a Cj-complex. Thus far, a definitive identifi¬ 
cation of the detailed local structure of these defects has 
yet to be made. 

2.2 Carrier lifetime studies Initial investigations of 
the minority earner lifetime (MCL) in 4H-SiC epitaxial 
materials were primarily carried out using optical tech¬ 
niques; time-resolved photoluminescence (TRPL) of the 
near-band edge emission at room-temperature, transient 
free carrier absorption (TFCA) and microwave photocon¬ 
ductivity decay (p-PCD). Using low-injection TRPL, 
Kordina et al. [38] found that the MCL exhibited a signifi¬ 
cant spatial variation across the wafer, with lifetime map¬ 
ping indicating variations from <0.5 ps to 2.1 ps, suggest¬ 
ing good material quality on portions of the wafer. These 
variations in lifetime were found correlated with the mor¬ 
phology of the wafer, suggesting that structural defects 
could affect the lifetime. Galeckas et al. [39] used TFCA to 
observe similar values for the carrier lifetime, albeit under 
higher excitation conditions. Evidence of contributions to 
the lifetime from surface and Auger recombination was 
also investigated. Spatial scanning provided TFCA lifetime 
maps that also reflected a correlation of structural defects 
with the lifetime. Using depth profiling, TFCA measure¬ 
ments [40] also provided the time-dependent depth distri¬ 
bution of the injected carrier concentration, thus enabling 
the determination of surface recombination velocities asso¬ 
ciated with the surface and with the nVn + -substrate inter¬ 
face. This demonstrated the extent to which the surface 
could play an important role in limiting the carrier lifetime, 
particularly in thinner layers. 

While the present focus is on point defects that limit 
the carrier lifetime, the observation of lifetime correlations 
with structural defects serves as a reminder that under the 
right conditions, other processes can also make significant 
contributions. TRPL mapping has shown, for example, that 
the lifetime can be influenced by structural defects in the 
substrate that propagate into the epilayers, and that correla¬ 
tion of the lifetime maps with X-ray topography studies as¬ 
sociates reduced carrier lifetimes with small angle grain 
boundaries containing threading edge dislocations [41, 42], 
However, while the observation that structural defects can 
affect the lifetime is well documented, their precise role is 
still unclear. For example, it has been reported [43] that the 
carrier lifetime in bulk wafers increases near regions of 
high defect density, due to the gettering of point defects by 
the structural defects. Similarly, the Z, concentration has 
been observed to decrease in epitaxial layers containing 
higher dislocation densities [44], Structural defects that in¬ 


volve local regions of a lower band gap SiC polytype, such 
as 3C inclusions in 4H, can also trap carriers. Stacking 
faults (SF’s), for example, made up of a few atomic layers 
of 3C within the 4H epilayers, effectively act as quantum 
wells [45] that can efficiently trap carriers and reduce the 
carrier lifetime in the region near the SF [46], As structural 
defects are localized, their influence on the bulk lifetime 
will depend in part on the fraction of the epilayer volume 
that they influence. Consequently, when there are high 
concentrations of structural defects, they can be effective 
in trapping carriers and can contribute significantly to lim¬ 
iting the measured lifetime. 

The first real attempt to combine DLTS with carrier 
lifetime measurements, in order to identify the point de¬ 
fects controlling the lifetime in 4H-SiC, was reported by 
Storasta et al. [47], who found a direct correlation between 
an increasing concentration of shallow boron centers and 
the decrease in the carrier lifetime measured by TRPL. Al¬ 
though this result has not been reproduced in the literature, 
the authors were careful to point out several specific diffi¬ 
culties with simply assigning boron as the lifetime killer 
based on these measurements. This also illustrates the fact 
that, while correlation is a highly valuable tool that aids in 
associating specific defects with observed material proper¬ 
ties, interpretations based on correlations alone must be 
made with caution. 

In addition to correlating DLTS signatures with growth 
and processing conditions, as described above, Zhang et al. 
[25] also found that variations in defect concentrations 
across the wafer were directly correlated with the carrier 
lifetime determined by TRPL. The results, shown in Fig. 6, 



Figure 6 Correlation between Z 1/2 and EH 6/7 concentrations and 
the carrier lifetime, for different sample positions along the 
gas flow direction. Reprinted with permission from Ref. [25], 
© 2003, American Institute of Physics. 
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indicated no obvious relationship between the carrier life¬ 
time and the concentration of shallow B, but did indicate a 
clear relationship between the decrease in the lifetime and 
the increase in the concentrations of both Z 1/2 and EH 6/7 . 
This result led to the conclusion that both defects were 
responsible for limiting the carrier lifetime. A similar cor¬ 
relation (measured by p-PCD and TRPL) between the life¬ 
time and both the Z 1/2 and EH 6/7 concentrations was ob¬ 
tained by Tsuchida et al. [48], who also noted a lack of 
correlation between the carrier lifetime and the D,-center 
or with the concentrations of several impurities. This 
behavior was also consistent with the later observations 
of Jenny et al. [49], who demonstrated that the high- 
temperature annealing of bulk, high purity 4H-SiC sub¬ 
strates to 2600 °C resulted in a dramatic increase in the 
carrier lifetime (<10ps to >3 ps) as well as a dramatic 
decrease in both the Z 1/2 and EH 6/7 defect concentrations. 
The above observations were all taken as evidence that 
both defects were controlling the carrier lifetime in the 
n-type epilayers. While this conclusion is quite reasonable 
given the experimental evidence, it is not the only possible 
interpretation. 

We have already seen that the concentrations of these 
two defects generally vary in the same way with changes 
in various growth or post-growth processing conditions, 
such as varying C/Si ratio, growth temperature, irradia¬ 
tion/implantation and annealing (see e.g. Figs. 3-6). This 
has been interpreted as reflecting the fact that both centers 
share the same intrinsic defect in their structure, most 
probably a V c . It follows then, that no matter what param¬ 
eter is varied to effect a change in the defect concentration 
(and thus in the carrier lifetime), the concentrations of 
these two defects will usually both respond in a similar 
way to this changing parameter. Even if only one of the 
two defects actually has any effect on the carrier lifetime, a 
correlation of lifetime vs. defect concentration will show a 
similar variation for both defects, simply because both de¬ 
fect concentrations always vary similarly in response to a 
changing growth or processing parameter. In effect - the 
lifetime is varying in response to a change in the concen¬ 
tration of one defect, the lifetime-killing defect, while the 
concentration of the second defect is simply varying with 
environment in the same way as the first defect, and has no 
effect on the lifetime. Consequently, it appears that both 
defects are correlated with the lifetime, when in fact only 
one is actually responsible for the reduction observed in 
the carrier lifetime. 

In Section 3.2 it is shown that this is, in fact, the case 
for 4H-SiC:EH 6/7 does not significantly affect the carrier 
lifetime, while Z 1/2 is identified as the dominant lifetime 
killer. While both of these defects have large electron cap¬ 
ture cross-sections and both are deep levels, it is shown 
that what distinguishes the two defects is that EH 6/7 has a 
small capture cross-section for minority holes. Conse¬ 
quently, Z 1/2 is an effective recombination center and mi¬ 
nority carrier trap while EH 6/7 is limited to the role of an ef¬ 
fective majority carrier trap. 


3 Identification of the lifetime killer 
3.1 Experimental Samples were grown by silane- 
based MOCVD on n + 4H-SiC substrates, with carrier con¬ 
centrations ~mid-10 14 cnT 3 and thicknesses ranging from 
9-104 pm [2], DLTS spectra were obtained in the 
usual way [2], while MCTS measurements employed opti¬ 
cally injected pulses from a 365 nm, light emitting diode 
instead of electrical injection. DLTS spectra were also 
measured on similarly grown 1.2 mm diameter p-i-n di¬ 
odes with n-layer thicknesses of 30-100 pm. Time- 
resolved photoluminescence measurements were carried 
out at 300 K and excited with either a frequency-tripled, 
Q-switched Nd:YAG laser (355 nm) or a frequency 
doubled, cavity-dumped, modelocked Ti: sapphire laser 
(355 nm, <0.1-1 MHz). The collected light was dispersed 
by a 0.22 m. double spectrometer and detected by a GaAs 
photomultiplier and time correlated single photon counting. 


3.2 Z 1/2 or EH 6/7 ? In order to investigate the proper¬ 
ties of these two defects, DLTS measurements were em¬ 
ployed to determine defect concentrations, while low- 
injection TRPL provided a measure of the minority earner 
lifetime. The defect concentrations are expected to be di¬ 
rectly related to the bulk MCL due to trapping at defects: 

4k = o’pVtf.Ak > (!) 

where er p is the minority carrier (hole) capture cross- 
section, v th the carrier thermal velocity, v th = (3 kT/m*) 1/2 , 
and N T the defect concentration. Carrier lifetime measure¬ 
ments, however, determine only an effective lifetime, r eff , 
which generally reflects contributions from capture and 
emission at defects, i.e. Shockley-Read-Hall recombina¬ 
tion (SRH) [50, 51], as well as from other competing re¬ 
combination processes, such as Auger recombination, sur¬ 
face recombination and radiative recombination: 


1 


Aff 



2"bulk (SRH) T Rad T Au ger f Sulf 


( 2 ) 


Because SiC is an indirect band gap material, the re¬ 
combination of free electrons with free holes is a very slow 
process. Consequently, the carrier lifetime is determined 
by the much faster competing processes. Hence, phonon- 
assisted radiative processes contribute only a small fraction 
to the total recombination rate. For low-to-moderate injec¬ 
tion conditions, Auger recombination is negligible. Thus, 
the bulk lifetime is dominated by SRH recombination. Sur¬ 
face recombination, however, cannot be ignored, particu¬ 
larly for thin epitaxial layers. Consequently, in order to re¬ 
late the bulk lifetime with the concentration of a particular 
defect using Eq. (1), the bulk lifetime must be extracted 
from the measured effective lifetime. This requires deter¬ 
mining the surface recombination rate, (r Surf )~', which can 
be reasonably approximated by [52] 


— = f4-+—T. 
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where d is the layer thickness, D the injection-dependent 
diffusion coefficient for 4H-SiC [53], and S represents an 
average surface recombination velocity for the layer. Con¬ 
sequently, one way to extract the bulk and surface lifetimes 
for a group of similar samples (i.e. assuming similar sur¬ 
face recombination) is to measure the dependence of the 
lifetime on layer thickness, d. A set of epitaxial layers was 
grown [2] for this purpose with a range of layer thick¬ 
nesses varying from 9 pm to 104 pm. 

DLTS and lifetime measurements were carried out on 
material taken from the center, middle and edge of each 
wafer. The variation in the concentrations of Z 1/2 , EH 6/7 and 
Ti centers are shown as a function of layer thickness for 
“middle” and “edge” samples in Fig. 7(a) and (b), respec¬ 
tively [54, 55], For the latter, the Z 1/2 and EH 6/7 concentra¬ 
tions show similar behavior, in agreement with the numer¬ 
ous reports discussed above. For the “middle” samples, 
however, there is a clear difference in the variations of the 
concentrations of these two defects with layer thickness. It 
can be noted that there is a general trend toward lower in¬ 
trinsic defect concentrations for thicker layers. While the 
cause of this is unclear, Fujihira et al. [56] have speculated 
that the stress between the epilayer and substrate leads to a 
higher defect concentration, and that trap formation may 
consequently be reduced in the thicker layers. It should 
also be noted that these variations in the defect concentra¬ 
tions for different thickness wafers do not conflict with the 
observed lack of any significant variation in defect concen¬ 
tration across the thickness of a single layer, recently re¬ 
ported by ul Flassan et al. [57]. In any case, the different 



Layer Thickness (pm) 

Figure 7 Variation of defect concentrations with layer thickness 
for samples obtained from the “middle” (a) and “edge” (b) re¬ 
gions of each wafer. (Adapted from Ref. [54].) 


dependences of the concentrations of these two defects on 
varying layer thickness, as in Fig. 7(a), provides an oppor¬ 
tunity to make a meaningful comparison between their 
measured defect concentrations and the carrier lifetime. 
The variation in concentration of Z 1/2 , EH 6/7 , shallow B and 
Ti centers are shown in Fig. 8(a) and are compared to the 
inverse of the effective carrier lifetime, as measured by 
low-injection TRPL. It is apparent from the figure that 
only the Z 1/2 concentration exhibits any clear correlation 
with l/r eff . This is, perhaps, a first indication that only one 
of these defects may limit the lifetime. 

The measured dependence of the inverse effective life¬ 
time on layer thickness may be fitted using Eqs. (1)—(3), 
with only SRH (bulk) and surface contributions. The result 
of this fit is shown as the open squares and dotted line in 
Fig. 8(b) and compared to the measured values (solid cir¬ 
cles). The fitted curve is not smooth because measured de¬ 
fect concentrations were used rather than requiring an ad¬ 
ditional fitting parameter. The best fit parameters were 
found to be S = 1250 cm/s and a p = 3.5 x 1(T 14 cm 2 [2], 
From this, bulk lifetimes could be extracted, and are plot¬ 
ted as a function of Z 1/2 concentration in the inset of the 
figure. The linear relationship again suggests that Z, /2 is the 
responsible defect. 



Layer Thickness (pm) 


o 

<o 

o 


Figure 8 (a) Variation of defect concentrations and the inverse of 
the TRPL decay time on layer thickness, and (b) fit (open 
squares) of the lifetime data (filled circles) using Eqs. (1) (3). 
The inset plots the extracted bulk lifetimes vs. Z 1/2 concentration. 
Reprinted with permission from Ref. [2], © 2006, American In¬ 
stitute of Physics. 
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Thus far, EH 6/7 has been ruled out as a lifetime-limiting 
defect based only on the lack of correlation between the 
defect concentration and the carrier lifetime within this 
particular group of “middle” samples. However, a less 
sample-dependent procedure to test the role of EH 6/7 in re¬ 
ducing the lifetime is preferable. This was carried out by 
DLTS measurements in a p-i-n diode structure [2, 54], 
Under a large forward bias, the drift region of the diode is 
flooded with a high concentration of minority holes. The 
filled electron traps Z 1/2 and EH 6/7 would be expected to 
trap these holes, changing their charge states and conse¬ 
quently reducing the strength of their DLTS signatures. 
Figure 9(a) shows a plot of DLTS signal strength vs. 
forward bias for both defects. Above 2 V the standard Z 1/2 
(-/+) signature is totally quenched and a (0/+) signal ap¬ 
pears - both resulting from the capture of minority holes 
by the doubly-occupied Z 1/2 . However, the EH 6/7 signal is 
not significantly affected by the large injected hole concen¬ 
tration, indicating that this defect must have a small cap¬ 
ture cross-section for minority holes. This observation has 
recently been confirmed by Reshanov et al. [58], as shown 
in Fig. 9(b). With a small hole capture cross section, EH 6/7 
cannot be responsible for trapping minority carriers and 
limiting the MCL. Therefore, Z 1/2 alone must be the lifetime 
killer. The carrier dynamics at this defect, which determine 
the MCL, will be considered in detail in a later section. 



Temperature (K) 

Figure 9 (a) Variation of the DLTS signals for Z 1/2 and EH 6/7 de¬ 
fects as a function of forward bias in a p-i-n diode. (Adapted 
from Ref. [54]). These observations were verified independently 
in (b) [58], © Trans-Tech Publications, 2009. 


3.3 Recent developments Recently, there has been 
further support of the view that Z 1/2 and EH 6/7 are related to 
the V c . Danno et al. [28] have determined that the forma¬ 
tion of these defects is independent of growth rate, and de¬ 
pends primarily upon C/Si ratio and growth temperature. 
The dependence upon the latter provided an activation en¬ 
ergy that was found consistent with a simple thermody¬ 
namic model for the carbon vacancy. The authors also ar¬ 
gued that the thermal stability of the defects was consistent 
with a V c , but not with C-interstitials. This view was sup¬ 
ported by Storasta et al. [59], who demonstrated that 
C-implantation and annealing led to a dramatic decrease in 
the Z 1/2 and EH 6/7 concentrations and a doubling of the car¬ 
rier lifetime, which was discussed in terms of the diffusion 
of implantation-induced C-interstitials and subsequent an¬ 
nihilation of carbon vacancies. It was later shown [60] that 
this approach resulted in thick epitaxial layers with a Z 1/2 
concentration <10 n cnT 3 , with resulting device structures 
that exhibited good conductivity modulation. By measur¬ 
ing the dependence of the carrier lifetime over a large 
range of Z, /2 concentrations, Danno et al. [61] have shown 
that below about 10 13 cnT 3 the carrier lifetime becomes 
sub-linear with Z 1/2 concentration, and at even lower con¬ 
centrations it becomes independent of Z 1/2 . At this stage, it 
is apparent that the contribution by this defect to the in¬ 
verse lifetime (e.g. Eqs. (1) and (2)) has become small 
compared to other sources, such as surface recombination 
or carrier capture by other defects. It should be noted, how¬ 
ever, that the concentration (10 13 cnf 3 ) at which this transi¬ 
tion occurs will be sample-dependent, as the surface condi¬ 
tion and the concentrations of other defects will clearly 
vary between different samples. While the defect or sur¬ 
face recombination that controls the lifetime in low-Z 1/2 
material may not compromise the lifetime sufficiently to 
be a major problem for lower-voltage bipolar devices, in 
high voltage switches, where the drift layer must be thick 
and the MCL must be very long, these remaining limits to 
the carrier lifetime can become problematic. 

As Z, /2 is the dominant lifetime killer when present in 
sufficient concentration, it is of interest to examine the 
processes associated with carrier capture at this defect, in 
order to gain a better understanding of its unique role and 
of the important parameters that determine the lifetime. 
This is considered in detail in the following section. 

4 Carrier dynamics at the Z 1/2 defect 

4.1 Introduction A study of the dominant processes 
controlling carrier capture and emission at the lifetime lim¬ 
iting defect in 4H-SiC epilayers provides an understanding 
of the key parameters that affect the lifetime, both in as- 
grown material and in device structures. The wide varia¬ 
tion in carrier lifetimes appearing in the literature can be 
understood to result from differences in the experimental 
conditions under which the measurements are carried out 
(e.g. injection level, temperature), differences in material 
quality (defect concentrations, surface condition), or even 
in the choice of experimental measurement technique, as 
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different measurements can be sensitive to qualitatively 
different lifetimes [62], Consequently, minority carrier 
lifetime measurements earned out at low injection in as- 
grown material can be very different from the lifetime in a 
device that has been subsequently fabricated from the same 
material and operating under high injection conditions and 
at elevated temperature. For carrier dynamics determined 
by SRH recombination at a given temperature, the lifetime 
of the minority carriers increases with injection level from 
the low-injection limit (r MCL ) dominated by minority car¬ 
rier capture, to the high injection limit (r HL ), where non¬ 
equilibrium electron and hole concentrations and lifetimes 
are essentially equal. Typical dependences of carrier life¬ 
time upon injection level, measured by TRPL, is shown in 
Fig. 10(a) and (b) for 160 pm and 52 pm thick n epilayers, 
respectively. A strong dependence on injection is observed 
between the low-injection and high-injection regimes. 
Measurements carried out by TFCA at the same position 
on the 160 pm thick wafer reflect about a factor of 2 longer 
lifetimes at high injection than the TRPL result, in agree¬ 
ment with theoretical comparisons between the two meas¬ 
urement techniques [62]. 

The effect on the carrier lifetime of external conditions, 
such as injection level and temperature, is approached ex- 



Figure 10 Injection dependence of the PL decay time for (a) 
sample A and (b) sample B. The PL decay time is compared in 
(a) to that determined by TFCA at the same location on the sam¬ 
ple. The solid lines are the results of fits of the data to carrier dy¬ 
namics simulations. Reprinted with permission from Ref. [62], 
© 2008, American Institute of Physics. 


perimentally and from a theoretical point of view, where we 
assume that the lifetime is determined by SRF1 recombina¬ 
tion at the Z m defect, so that Auger and surface recombina¬ 
tion are considered negligible. The experiments were ar¬ 
ranged to approximate these conditions: Surface recombina¬ 
tion is minimized by limiting the measurements to relatively 
thick layers, as in Fig. 10, and injection levels are limited to 
<10 18 cnT 3 , where Auger recombination is negligible [40], 

The capture and emission of carriers at defect sites is 
usually described theoretically by the SRH recombination 
model [50, 51]. However, this picture assumes a simple de¬ 
fect structure (e.g. two possible charge states - fdled and 
empty) and steady state excitation conditions, and is not di¬ 
rectly applicable to a complex defect such as Z 1/2 and tran¬ 
sient excitation. Although much work has been done over 
the years to extend the SRH model to multi-level systems 
and to transient conditions [63-69], the problem is not 
soluble in closed form. Consequently, much of this work 
has been limited to very high or very low injection, where 
the problem has relatively simple solutions. This is not a 
useful approach in the present case, where the variation in 
the lifetime for arbitrary injection level and temperature is 
of interest, and where the defect has a more complicated 
internal level structure. Consequently, we will take the ap¬ 
proach of employing numerical carrier dynamics simula¬ 
tions to explore the effects of injection level and tempera¬ 
ture on the carrier lifetime and to compare the calculated 
results to experimental measurements. 

The model employed for the simulation is based on the 
fact (see e.g., Section 2) that Z 1/2 has been shown to be a 
-U defect that can trap two electrons, leading to three pos¬ 
sible charge states; doubly occupied (“-”), singly occupied 
(“0”) and empty (“+”). Because of the large lattice relaxa¬ 
tion around the doubly-occupied Z state (see e.g. inset in 
Fig. 2), the additional repulsive energy associated with the 
two-electron state is more than offset by the reduction in 
lattice energy. Consequently, the second electron to be 
captured is bound more tightly than the first, and the two- 
electron Z state is the equilibrium charge state in n-type 
material [3], Transitions between Z and Z° (—/0), or be¬ 
tween Z° and Z + (0/+), occur through the capture and emis¬ 
sion of electrons and holes, and thus affect the carrier life¬ 
time directly. These charge-state transitions are shown 
schematically in Fig. 11. The energies £j and E 2 in the fig¬ 
ure represent thermal activation energies for the (—/0) and 
(0/+) transitions, respectively. As a result of the large lat¬ 
tice relaxation, the ( /0) thermal emission barrier also in¬ 
cludes a capture barrier, E b , as shown in the inset of Fig. 2. 
The R,j and R pj (j= 1, 2) in the figure represent net electron 
and hole capture rates for the two transitions 1 and 2, cor¬ 
responding to ( /0) and (0/+), respectively. Solid lines in 
the figure represent carrier capture, while dotted lines rep¬ 
resent the thermal emission of carriers. By employing this 
model to represent the Z 1/2 defect in carrier dynamics simu¬ 
lations, the dependence of the carrier lifetime for arbitrary 
variations in injection level and temperature can be inves¬ 
tigated and compared to experiment. 
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Figure 11 Schematic level scheme for carrier capture and emis¬ 
sion at the Z 1/2 defect. The three charge states of Z 1/2 (-, 0, +) are 
indicated, resulting in two transitions: 1 for (-/0) and 2 for (0/+). 
Transitions due to carrier capture or emission are indicated by 
solid or dashed arrows, respectively, with corresponding trap 
depths £j and E 2 . The net carrier capture rate for each transition is 
indicated by the R nj and R pj . Reprinted with pennission from Ref. 
[62], © 2008, American Institute of Physics. 

4.2 Calculation parameters For the simulation to 
reflect the experimental situation with reasonable accuracy, 
the model employed must be a good approximation to the 
experimental conditions, and an accurate set of material 
parameters must be available. There is ample evidence that 
Z 1/2 dominates the carrier lifetime, as long as the defect 
concentration is at a non-negligible level [61]. Conse¬ 
quently, SRH recombination at this defect appears to be a 
reasonable picture, as long as Auger and surface recombi¬ 
nation can be neglected. The extent to which the calculated 
results agree with experimental measurements is sensitive 
to the particular set of material parameters employed in the 
calculation. There are some general trends, however, that 
are noteworthy: The variation of the lifetime with injection 
level tends to be most sensitive to carrier capture cross- 
sections and to defect concentrations, while the tempera¬ 
ture dependence is most strongly affected by the thermal 
ionization energies E x and E 2 , as these enter exponentially 
in determining the emission rates of the carriers from the 
deep centers. 

Many of the parameters required in the simulations 
have been measured and are available in the literature. Ef¬ 
fective masses for 4H-SiC have been determined by Son 
et al. [70], while trap depths and electron capture cross- 
sections for Z 1/2 were reported by Hemmingsson et al. [3]. 
Since the material parameters for Z, and Z 2 are similar, Z 1/2 
is treated here as a single defect, with averages of the simi¬ 
lar parameters for Z, and Z 2 reported in Ref. [3] used in the 
calculation. While Hemmingsson et al. have measured the 
temperature dependent (0/-) electron capture cross section 
cr„|, the (+/0) electron capture cross section a„ 2 reported 
by this group and by several others falls into a wide range: 


Table 2 Parameters used in numerical solution of Eqs. (5)—(7). 


parameter 

symb value 

Ref. 

-/0 trap depth” 

E l 

0.62 eV 

[3] 

0/+ trap depth” 

e 2 

0.47 eV 

[3] 

0/- capture barrier” 

E b 

0.07 eV 

[3] 

0/- e~ capt.cross section” 

°n\ 

1.5 x 10 15 e -°- 07/kT cm 2 

[3] 

+/0 e~ capt. cross section 

<*nl 

*(2-4) x 10 15 cm 2 

c 

-/0 h + capt. cross section 

°pl 

3.5 x 10~ 14 cm 2 

[2] 

0/+ h + capt. cross section 

°p2 

*(1-2) x 10~ 14 cm 2 

c 

Zj/Z 2 concentration 13 

N 

10 11 — 10 14 cnT 3 


equilibrium carrier conc. b 

«o 

10 14 -10 15 cnT 3 



3 Average of Zj andZ 2 values determined in Ref. [3]. 
b Typical range of concentrations found in n epilayers. 
Fitted value from simulations. 


<t „2 « (1-20) x 10- 15 cm 2 . This uncertainty results from the 
fact that these values were obtained as an extrapolation of 
measured data, and whether or not a temperature depend¬ 
ence had been assumed for the capture cross section. Be¬ 
cause of this uncertainty, the value of a„ 2 is estimated here 
by varying a„ 2 to fit the experiment. The (-/0) hole capture 
cross section a pl has recently been measured [2], but no 
measurement is available for the (0/+) hole capture cross 
section, a p2 . Consequently, this parameter is also treated as 
a variable to be fitted to the experimental results. The other 
significant parameters in the calculation are the defect con¬ 
centration N and the equilibrium carrier concentration n 0 , 
which are sample-dependent. In determining a n2 and a p2 , 
measured values are used for n 0 . However, since DLTS 
measurements were not available for these samples, the de¬ 
fect concentrations are determined from the measured low- 
injection minority carrier lifetime: r MCL = (a pX N) Some of 
the important parameters employed in the simulations are 
summarized in Table 2. With reasonable fitted values for 
these two capture cross-sections, generalized simulations 
can be carried out over a large range of injection level or 
temperature. 

4.3 Experimental TRPL measurements were de¬ 
scribed in Section 3.2. A flow-through optical cryostat or a 
temperature-regulated heat stage was used for temperature- 
dependent measurements. Transient free carrier absorption 
measurements employed the Nd:YAG laser (Section 3.2) 
as a pump source, a CW, 1310 nm diode laser (0-5 mW) 
probe beam and a 0.3 mm diameter InGaAs detector 
matched to a transimpedance amplifier. Three n epilayers 
(ra 0 ~ mid 10 14 cm 3 - low-10 15 cm -3 ) were investigated: 
Samples A (3"-wafer, 160 pm thick) and B (*lcm 2 , 
52 pm thick) were used for injection dependence studies, 
and sample C («1 cm 2 , 160 pm thick) for temperature de¬ 
pendence measurements [62]. 

4.4 Carrier dynamics simulation Referring to 
Fig. 11, the SRH carrier lifetime is determined by the cap¬ 
ture and emission of electrons and holes. In general, if 
n and p are electron and hole concentrations defined by 
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n(t) = n 0 + bn(t); p(t) = p 0 + bp(t), where the “0” subscripts 
refer to equilibrium values and 8w, S p represent their 
variations from equilibrium, and with p 0 * 0 in n-type ma¬ 
terial, instantaneous carrier lifetimes can be defined 
as T n (t) = [-(AbnlAt)lbn)) ', v p (t) = [—(d5/?/dr)/5/?] 1 . Since 
the time-varying PL intensity is then of the form 
l„(t) = y[n 0 + bn(t)\ 8 p(t), the decay time measured by 
TRPL can be written r PL (t) = [-(d/p L (t)/d?)//pL(t)] _l . In terms 
of the instantaneous electron and hole lifetimes, this can be 
shown to take the form 


1 1 1 6 n(t) 

7Vl(0 r p (t) r„(t)[n 0 +bn(t)_ ' 


(4) 


At very low injection (bn —> 0), the PL decay time is 
identical with the hole, or minority carrier lifetime, r MCL , 
while at very high injection (bn > n 0 and r„ « T p = r HL ), 
r PL —» r HL /2 due to the bimolecular recombination kinetics 
at high injection [62], Similarly, the time-varying concen¬ 
trations of the three Z, /2 charge states, Z k (t), k = (+, 0, -,) 
can be written as Z*(?) = + SZ*(7). The 8Z*(t), bn(t) and 

8 p(t) constitute the five time-dependent variables defining 
the problem, of which two can be eliminated, since the to¬ 
tal defect concentration N is fixed and charge is conserved: 

^bZ k = 0, bn -bp = 8Z° + 2SZ + . (5) 

k 


Eliminating two of the defect concentrations, the time 
dependence of the carrier and defect concentrations are 
then determined by 


d8n 
d t 
dbp 
d t 


~(R„\ +^» 2 )> 


-( R p\ +R pi)’ 


d8Z~ 
d 1 


= R„i-R pl 


( 6 ) 


where each of the net capture rates (e.g. Fig. 11) are de¬ 
termined from the difference between capture and emission 
for each transition, and take the form [62]: 

R„\ = a„ i {(« 0 + 8 n) [8p -bn-2bZ~ + Z 0 °] 
-n x [Z-+bZ~}}, 

R„i =a n2 {{n 0 +bn)[bn-bp + bZ~ +Zl] 

- n 2 [bp -bn-2bZ~ +Z 0 0 ]}, 

R p i = ° P \[(Po + 8 p) (Z~ + bZ ~)] 

-p, (bp -bn- 2 bZ~ + Z 0 °), 

Rp 2 = a p2 [(Po + 8p) (8 p - bn - 2bZ~ + Z 0 °) 

—pi (bn - bp + bZ~ + Z 0 + )]. 


Here, the a nj and a pj are electron and hole capture coef¬ 
ficients; e.g. a nJ = ct^v^h, with electron capture cross sec¬ 
tion ct„j and and hole capture cross section a pj , respectively, 


and v th np = (3 kT/m* h ) V2 are the average carrier thermal ve¬ 
locities. The quantities and p ; are given by [50] 
nj = A c exp [-E/kT] and pj = N v exp [-(E s -Ej)/kT], with 
N C (T) and N V (T) the effective density of states for the con¬ 
duction and valence bands, respectively and the £] denote 
the energy separations from the conduction band of the two 
defect levels (-/0) and (0/+). 

The primary independent variables of interest are the 
injection level and the temperature. For a short excitation 
pulse, the injection level g enters the calculation as an ini¬ 
tial condition for the nonequilibrium carrier concentra¬ 
tions: 8«(0) = 8/2(0) = g. The temperature variations of the 
nonequilibrium carrier concentrations and their lifetimes 
are determined by the temperature dependent parameters in 
the calculation: The n. and pj vary as 7’ 3/2 exp \-EJkT\ and 
dominate over the temperature dependence of the a ly and 
a pj , which vary as T' 11 , in addition to any temperature 
variation of the capture cross sections. For T> 100 K, the 
equilibrium carrier concentration is in the saturation range 
[71], so the equilibrium carrier concentration w 0 may be 
taken as temperature-independent. 

4.5 Simulation results The carrier dynamics simula¬ 
tions provide the time dependence of the nonequilibrium 
carrier concentrations and the three defect charge state 
concentrations after a short, above-gap optical pulse, 
and for varying injection level and temperature. In order 
to probe the general features of carrier recombination at 
Z 1/2 , it is first necessary to fit experimental results to 
Eqs. (4)—(7) in order to obtain consistent values for the 
two capture cross-sections <j„ 2 and <j p2 , so that a full set of 
parameters (Table 2) is available. The resulting fit of the 
dependence of TRPL decay time on injection level is 
shown for samples A and B as the solid lines in Fig. 10(a) 
and (b). The fitted values for the two capture cross sections 
for both of these measurements, as well as for the fit of the 
temperature dependence of the TRPL decay time for sam¬ 
ple C, discussed later in this section, all fell within the 
range a„ 2 ~ (2-4) x 10~ 15 cm 2 and <j p2 ~ (1-2) x 10~ 14 cm 2 . 
With a consistent set of parameters, the general features of 
carrier dynamics after pulsed excitation can be considered. 
While the results are dependent on the material parameters 
employed, e.g. defect concentration and doping level, there 
are several qualitative features apparent in Table 2 that will 
have a dominating effect on the results. First, the capture 
cross section cr nl is much smaller («1(L 16 cm 2 at room tem¬ 
perature) than the other three cross-sections, so that (-/0) 
electron capture or emission will be slower, and can be a 
rate-limiting process. Also, the large depth of the defect 
levels constrains electron emission to be negligible 
(«!, n 2 bn, bp) for T< 500 K. Finally, the large band gap 
of 4H-SiC results in very high barriers (>2.5 eV) for hole 
emission from these traps. Consequently, hole emission is 
insignificant (p\,p 2 ~ 0) for temperatures <1300 K. 

Using the parameters in Table 2, general trends in the 
carrier dynamics can now be examined. Initially, at room 
temperature the defects are essentially all in the Zr state. 
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Figure 12(a) shows the calculated time dependence, after a 
short excitation pulse, of the nonequilibrium carrier con¬ 
centrations, 5 n(t) and S p(t), and the three defect charge 
state concentrations Z k (t), under high injection conditions, 
g = 1 x 10 16 cnr 3 , and with typical defect concentration 
N = 1 x 10 13 cnT 3 and equilibrium carrier concentration 
n 0 = 4 x 10 14 cm -3 . The results suggest that the return of the 
system back to equilibrium follows four main stages of de¬ 
cay, which are indicated by the Roman numerals in the 
figure. Stage 1 involves rapid initial hole capture by the 
equilibrium Z state, on a time scale ( a pX gY 1 « 240 ps, 
which rapidly depletes the Z concentration and builds up 
the Z° and Z 3 populations. 

Stage 11 is a “high-injection stage”, where 8n, 5 p > n 0 , 
N, and therefore 5 n ~ 8p, so that the electron and hole 


lifetimes are approximately equal: r„ « r p = t H l- This phase 
of the carrier recombination can be understood by first 
considering the behavior of the charge state populations for 
the case where T < 500 K, so that thermal emission is neg¬ 
ligible (although the result is easily extended to higher 
temperature [62]). The time dependence of the Zr and Z 3 
concentrations are then determined by: dZ7d? = R nl -R pl 
and dZVdt = R p2 -R„ 2 , which can be written, using Eq. (7) 
in the limit of high injection, as: 

—— + a ,5« Z” ( t ) = a nl 8n Z° (t ) , 

d ' ( 8 ) 

d7+ 

- \-a 2 8nZ + = a 2 8nZ° . 

At " 2 p 2 



0 5 10 15 

Time (psec) 


Figure 12 Simulated time dependence of the carrier concentra¬ 
tions and defect charge states following a short, above-gap optical 
pulse, for (a) 1 x 10 16 cnT 3 initial injection and (b) 1 x 10 14 cnT 3 
initial injection. In (a) four stages of decay of the carrier system 
are indicated, as discussed in the text. The corresponding evolu¬ 
tion of the PL and carrier lifetimes are shown in (c). Reprinted 
with permission from Ref. [62], © 2008, American Institute of 
Physics. 


The characteristic times over which Z~ and Z 3 (and thus 
Z°) vary, (a pX 8ny' and (a„ 2 8n)\ respectively, are short 
at high injection: 240 ps and 1.5 ns, respectively, at 
1 x 10 l6 cm 3 injection. For much longer times, the Z\t) 
approach their constant, quasi-steady-state (QSS) values 
(/ —> oo, dZVd/ —> 0). Therefore, the Z^{f) concentrations 
remain slow-varying at their QSS levels as long as the de¬ 
caying carrier concentration is still in the high-injection re 
gime. This can be seen in the slow-varying charge state 
concentrations in Fig. 12(a) within Stage II. Under QSS 
conditions, Eq. (8) becomes an algebraic equation, which 
can be solved (with ZZ k = N) for the individual QSS 
charge state concentrations in terms of N and the capture 
coefficients. 

As the decaying earner concentration approaches « 0 , 
the high injection limit no longer applies, the electron and 
hole lifetimes begin to diverge and the decay enters into 
Stage III, a transition between high and low injection. Here 
the depleted Z population begins to increase via 0/- elec¬ 
tron capture, while Z 3 decreases via +/0 electron capture. 
The increasing Z concentration leads to an acceleration of 
—/0 hole capture, resulting in a rapid decrease in the hole 
lifetime. Stage III is essentially identical to the decay at 
low initial injection, g= 1 x 10 14 cnT 3 , shown in Fig. 10(b), 
except for a rigid time shift of several microseconds, which 
roughly corresponds to the duration of Stage II under high 
injection. 

The last stage of decay, Stage IV, is characterized by 
low-injection. Here, 8p <§ 8n <§ n 0 and Z ~ N. The concen¬ 
tration of the few remaining nonequilibrium electrons is 
identical to the concentration of Z°, and decays in the low- 
injection limit with lifetime r„ —> ( a nl n o y\ while the hole 
lifetime is the minority carrier lifetime, r p —> (a pl N) '. 

The carrier dynamics in Fig. 12(a) can also be de¬ 
scribed by the time dependence of the instantaneous carrier 
lifetimes, as shown in Fig. 12(c). Here, Stage II (r„» T p ) 
and Stage IV (r„, T p « constant) are quite apparent. The PL 
decay time is also shown, with high-injection behavior 
(tpL ~ r HL /2) at early times and low-injection behavior 
(r PL « Tp) at later times, after the nonequilibrium earner 
concentration has decayed substantially. 
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From Fig. 11 it is clear that we can write 
d5«/dt = —[a„ t nZ° + a n2 nZ']. Under high injection condi¬ 
tions, where Z k » Zq SS and n ~ 5 n, the carrier lifetime 
can be written as r„' = r H ' L = -[(d5n/d/)/8n] = —[a Bl ZQ SS + 
a„ 2 Zjj SS ], where the Z^ ss are determined in terms of N and 
the capture coefficients, as noted above. Consequently, 
with r MCL = ( a pt N)', the ratio of the high and low injection 
lifetimes, (r HL /r MCL ) can be determined in terms of the cap¬ 
ture coefficients [62]: 


(Thi/^MCl) — 


1 + _p! 


( 9 ) 


For a„, <g a nl , a pl , a p2 , appropriate for Z 1/2 , this simplifies to 
%/tMCL ~ a P \{a „2 + a p2 )/a„ 2 a p2 . Using the parameters in Ta¬ 
ble 2, we find ( r HL /r MCL ) a 10. Incorporating the additional 
factor of 1/2 for the high-injection PL lifetime (r PL a r HL /2), 
the same ratio for PL decay times becomes (r HL /r MCL ) PL a 5, 
which agrees quite well with the measured high- and low- 
injection PL lifetimes in Fig. 10. 

Simulations were also carried out assuming that both 
traps, Z 1/2 and EH 6/7 , contribute significantly to the recom¬ 
bination. While the electron capture cross section for EH 6/7 
has been measured as cr„ 3 = 5 x 1(L 15 cm 2 at room tempera¬ 
ture [72], there is no measurement available for the hole 
capture cross section, a p3 . Since DLTS measurements in 
forward biased p—i—n diodes indicate this parameter to be 
small [2, 58], a value of 1 x 1CT 17 cm 2 was adopted for the 
simulations. The results, similar to Fig. 12(a), indicated lit¬ 
tle effect on the carrier lifetimes until the EFl 6/7 concentra¬ 
tion was increased far beyond what is observed in 
as-grown material. It is apparent that the small hole capture 
cross-section for this defect limits its role as an effec¬ 
tive recombination center, although the large electron cap¬ 
ture cross section and large thermal ionization energy 
(~ 1.6 eV) make it an effective trapping center for majority 
carriers. 

The change in the carrier lifetime with varying tem¬ 
perature was also considered, both experimentally and 
through carrier dynamics simulations. The calculated car¬ 
rier lifetime is affected by temperature in two ways; 
through the temperature dependence of the capture and 
emission rates and through the temperature dependence of 
the equilibrium distribution of the three Z 1/2 charge states. 
The latter is taken into account in the simulation as a tem¬ 
perature dependent initial condition in the solution of 
Eq. (6). It was noted earlier that the temperature depend¬ 
ence of the carrier lifetime controlled by SRH recombina¬ 
tion at the Z, /2 defect will be dominated by the exponential 
dependence associated with thermal emission from deep 
traps, and that hole emission will be negligible except at 
the highest temperatures. Consequently, electron emission 
is expected to dominate the temperature dependent lifetime. 
Since the thermal emission of trapped carriers tends to re¬ 
populate the decaying free carrier concentration, enhanced 


carrier emission due to increased temperature is expected 
to lead to a corresponding increase in the carrier lifetime. 

This increase is also observed experimentally [38, 62], 
In Fig. 13(a), the temperature dependence of the PL decay 
time, measured for sample C (open squares), is compared 
to the results of the simulation (filled circles). The calcula¬ 
tion reproduces the observed rapid increase in lifetime 
above 500 K reasonably well. The resulting fitting param¬ 
eters were cr„ 2 = 3 x 1CL 15 cm 2 and a p2 = 1.4 x 1(L 14 cm 2 , 
which are consistent with the corresponding parameters de¬ 
termined previously by fitting the injection dependences of 
the lifetime for samples A and B. The experimental points 
are compared in the figure to data reported by Kordina et al. 
[38], and agree reasonably well over their common range 
of temperatures. However, the measurements do not reflect 
the sharpness of the “knee” in the simulation near 500 K. 
This may well result from limitations in the model of 
Fig. 11, which assumes contribution from only a single de¬ 
fect. The more gradual rise in the lifetime with temperature 
below 500 K is consistent with a small contribution to 
the lifetime from other defects with lower thermal activa¬ 
tion energies than Z 1/2 , such as the boron hole trap 
(E v + 0.27 eV), which is filled in n-type material and avail¬ 
able to trap nonequilibrium minority holes. 
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Figure 13 (a) Measured temperature dependence of the PL life¬ 
time for sample C (open squares) compared to the simulated 
behavior (solid circles). Measurements reported in Ref. [38] 
are shown as filled diamonds, (b) The simulated temperature de¬ 
pendence of the three Zj/Z 2 charge state concentrations for 
t ~ 0.5 ps. Reprinted with permission from Ref. [62], © 2008, 
American Institute of Physics. 
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The dominant temperature dependence responsible for 
the rapid rise in lifetime above 500 K was found to result 
from electron emission from Z°, which is expected, as this 
transition involves the smallest of the two activation ener¬ 
gies. The nature of the transition can be seen clearly from 
the simulation: The simulated temperature dependence of 
each of the three Z, /2 charge states is shown in Fig. 13(b), 
for i=0.5 ps after the exciting pulse. It is apparent from 
the figure that as temperature increases toward 500 K, the 
concentration begins to increase, while the Z° concentra¬ 
tion decreases exponentially, with an activation energy 
near that of E 2 . These are the hallmarks of a 0/+ transition 
resulting from the thermal emission of electrons from Z°. 
The much smaller Zr concentration also decreases, with 
thermal activation energy £j, reflecting -/0 electron emis¬ 
sion. As temperature is increased further, Z° and Zr concen¬ 
trations continue to decline until, at very high temperatures, 
hole emission Z° —> Zr (> 1300 K) and ZZ —f Z° (>1600 K) 
become significant, and these two charge states begin to 
increase in population with increasing temperature. 

5 Summary In the preceding, we have followed the 
progress toward understanding the identity, the local struc¬ 
ture and the recombination dynamics of the defects that 
control the carrier lifetime in n-type 4H-SiC epitaxial lay¬ 
ers. Over the course of about a decade, the results of many 
measurements concluded that two defects, Z 1/2 and EH 6/7 , 
both dominant electron traps in DLTS spectra, were re¬ 
sponsible for controlling the carrier lifetime. Over about 
the same time period, a progression of DLTS studies were 
carried out with varying growth temperature, C/Si ratio, 
thermal anneal temperature, electron irradiation and ion 
implantation, concluding that the structure of both of these 
defects involved a carbon vacancy. 

However, later studies showed that the observed corre¬ 
lations between the carrier lifetime and both Z 1/2 and EH 6/7 
could result from the fact that the two defects, sharing a 
similar structure, tend to vary in concentration in a similar 
manner in response to changes in growth or post-growth 
processing parameters. Thus, the same correlation of both 
defects with the lifetime would be observed, even if only 
one of the two defects had a significant impact on the car¬ 
rier lifetime. This possibility was confirmed experimen¬ 
tally: Individual cases were observed where only Z 1/2 was 
correlated to the carrier lifetime, and DLTS studies of for¬ 
ward biased p-i-n diode structures indicated that EH 6/7 
has a small cross capture section for minority holes, and 
cannot be the dominant minority carrier trap. Consequently, 
it was concluded that Z 1/2 alone is the lifetime killer in this 
material. 

This has led to some interest in the details of the re¬ 
combination processes occurring at the Z 1/2 center. To this 
end, transient carrier decay measurements have been car¬ 
ried out over a wide range of carrier injection levels and 
temperatures, and were compared to the behavior calcu¬ 
lated for this complex defect. Since the carrier dynamics 
problem has no closed form solution for arbitrary injection 


level and temperature, simulations were performed using 
material parameters available in the literature. Two capture 
cross-sections that have not been determined uniquely 
were employed as fitting parameters, and the resulting val¬ 
ues were found consistent for three independent sets of 
lifetime measurements. Using a consistent set of param¬ 
eters, the generalized simulations provided insight into the 
details of carrier recombination at Z 1/2 . The response of the 
system to a short, above-gap optical pulse was observed to 
follow four distinct stages of decay for the nonequilibrium 
electron and hole concentrations and the populations of the 
three Z I/2 charge states. A high-injection phase, following 
the rapid capture of minority holes by the equilibrium Z 
state, was characterized by slow-varying, quasi-steady- 
state defect charge state concentrations, as both earner 
types decayed at the high-injection limit, r HL . A transition 
to the low-injection regime followed, as the decaying non¬ 
equilibrium carrier concentration approached the doping 
level. The temperature dependence of the carrier lifetime 
was found to be relatively weak below 500 K, but in¬ 
creased rapidly above this temperature, as electrons 
trapped in the Z° charge state are thermally emitted into the 
conduction band. At higher temperatures, contributions 
due to electron emission from the deeper Zr state occur. 
Due to the very large emission barrier, hole emission 
does not become significant until very high temperature 
(>1300 °C). 
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